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Abstract. This paper provides a rigorous and delicate analysis for exponential decay of 
Jacobi polynomial expansions of analytic functions associated with the Bernstein ellipse. 
Using an argument that can recover the best estimate for the Chebyshev expansion, we 
derive various new and sharp bounds of the expansion coefficients, which are featured 
with explicit dependence of all related parameters and valid for degree n > 1. We demon- 
strate the sharpness of the estimates by comparing with existing ones, in particular, the 
very recent results in 1381 SIAM J. Numer. Anal., 2012]. We also extend this argument 
to estimate the Gegenbaucr-Gauss quadrature remainder of analytic functions, which 
leads to some new tight bounds for quadrature errors. 



1. Introduction 

The spectral method employs global orthogonal polynomials or Fourier complex expo- 
nentials as basis functions, so it enjoys high-order accuracy (with only a few basis functions), 
if the underlying function is smooth (and periodic in the Fourier case). The convergence 
rate 0{n~ r ), where n is the number of basis functions involved in a spectral expansion and 
r is related to the Sobolev-regularity of the underlying function, is typically documented 
in various monographs on spectral methods [18l HH [14l [H [2H [35J LZJ [3 [24l [32]. It is also 
widely appreciated that if the function under consideration is analytic, the convergence 
rate is of exponential order 0(q n ) (for constant < q < 1). However, there appears very 
limited discussions of such error bounds (mostly mentioned, but not proved) in [Mj [35] [7] . 
Indeed, as commented by Hale and Trefethen [23] , the general idea of such convergence goes 
back to Bernstein in early nineties, but such results do not appear in many textbooks or 
monographs, and there is not much uniformity in the constants in the upper bounds. 

An important result in Bernstein [5] (1912) (also see [28]) states that u is analytic on 
[—1,1], if and only if 

sup lim yE N (u) = -, E N (u)= inf ||w-w||oo, 

JV->oo p vEPn 

where Pjy is the polynomial space of degree no more than N, and p > 1 is the sum of the 
semi-axes of the maximum ellipse £ p with foci ±1, known as the Bernstein ellipse, on and 
within which u can be analytically extended to. One immediate implication is that the best 
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polynomial approximation in the maximum norm enjoys exponential convergence. A more 
precise estimate for the Chebyshev expansion can be found in various approximation theory 
texts (see e.g., [HJ Theorem 3.8] and Theorem 5.16]): 

2M M 
\u%\<— , Vn>0; h-S^K {p _ 1)pN , (1.1) 

where M = max 2g £ {m„ } are Chebyshev expansion coefficients of u, and S^u is the 

partial sum involving the first N + 1 terms. One also refers to [33] [T2j [31] [6] [29] [36] [37] 
and the references therein for verification/description of exponential convergence of Fourier, 
Chebyshev or Legendre expansions. We remark that Gottlieb and Shu et al [20] [19] studied 
exponential convergence of Gegenbauer expansions (when the parameter grows linearly with 
n) in the context of defeating the Gibbs phenomenon. 

Here, we particularly highlight that a very recent paper of Xiang [38] provided a simple 
approach to obtain the bounds for Jacobi expansion coefficients of analytic functions on and 
within the Bernstein ellipse £ p : 



K' P \ < tW 2 ^ where = —b \ u{x)J^{x)u^{x)dx. (1.2) 

P (P - L ) V In Jn J -I 

Here, { J"'^}(a, (3 > —1) are Jacobi polynomials mutually orthogonal with the weight func- 
tion uj a ' f3 (x) = (1 - x) a (l + x)P and with the normalization factor ^ (cf. ([2T5])). The key 
step is to insert the Chebyshev expansion u{x) = ^2JLo ^ L jTj{x) into the Jacobi expansion 
coefficients and rewrite 



u n ' — 

7" j=o 



, oo -i 

jjE^ / T j (x)J^(x)u,^(x)dx, 

1. „-_n J — 1 



so the bound for the Chebyshev coefficient in (|1.1[) could be used. 

The first purpose of the paper is to take a different approach to derive sharp estimates for 
general Jacobi expansion of analytic functions. The assertion of sharpness is in the following 
sense: 

(i) The bound for general Jacobi case is tighter than (|1.2[) (see Remark |2.3|) . 

(ii) Refined estimates can be obtained for Gegenbauer expansion (a = /? > — 1), 
Chebyshev- type expansion (a = k — 1/2,(3 = 1 — 1/2 for non- negative integers 
fc, I), and Legendre-type expansion (a = k, f3 — I for non-negative integers k, I). The 
argument can recover the bounds known to be the sharpest (e.g., the Chebyshev 
case), and some obtained estimates are new and significantly improve the existing 
ones (see e.g., Remark [23]) . 

A second purpose of this work is to extend the argument to analyze Gegenbauer-Gauss 
quadrature of analytical functions. Recall that the remainder of Gauss-quadrature with the 
nodes and weights {xj, Uj} n =1 , takes the form (see e.g., [13]): 

u(x)uj(x) dx — ^^u(x.j)ujj = — (p " u{z)dz, (1.3) 

-1 j ^ £p *n\) 

where are the zeros of p n (x), orthogonal with respect to the weight function u>(x), 

and 

qn ( Z ) = U lp -^^dx. (1.4) 

Z ./ i Z X 
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The estimate of quadrature errors has attracted much attention (see e.g., JTJJ [Til] [3J [T7J 
IT5I [TBI |2"51 US]). Among these results, intensive discussions have been centered around the 
Chebyshev case and its family, e.g., Chebyshev of the second kind, but with very limited 
results even for Legendre-Gauss quadrature (see e.g., [HI HI]). In fact, the analysis heavily 
relies on the availability of explicit expression of p n {z) on £ p . Armed with a delicate estimate 
of q n (z) (in the first part of the paper) and the explicit formula of Gegenbauer polynomial 
in [39], we are able to derive sharp bound for the Gegenbauer- Gauss quadrature errors. 

We remark that there has been much interest in estimating spectral differentiation errors 
of analytic functions. Tadmor [34] first attempted to estimate the aliasing errors to verify 
exponential convergence of Fourier and Chebyshev spectral differentiation with a different 
assumption on analyticity. The results for analyticity characterized by the Bernstein ellipse 
include Reddy and Weideman [30] for Chebyshev case, and Xie, Wang and Zhao [39] for 
Gegenbauer spectral differentiation. It is also interesting to point out that Zhang [40} [41] [42] 
studied superconvergence of spectral interpolation and differentiation. We stress that the 
analysis apparatuses and arguments in this pipeline are different from these in this work. 

The rest of this paper is organized as follows. In Section [21 we provide sharp bounds 
for general Jacobi expansions of analytic functions, followed by some refined results for 
Chebyshev-type and Legendre-type expansions. In Section [3] we extend the argument to 
analyze Gegenbauer-Gauss quadrature errors. In the final section, we provide results to 
show the sharpness of the bounds by comparing them with existing ones. 

2. Sharp bounds for Jacobi expansions 

We derive in this section sharp bounds for Jacobi expansions of functions analytic on and 
within the Bernstein ellipse £ p . 

2.1. Preliminaries. It is known (see e.g., [12]) that the Bernstein ellipse is transformed 
from the circle 

C p = {w = pe w :9e [0,2tt]}, p > 1, (2.1) 
via the conformal mapping: z = (w + w )/2, namely, 

£ p := {ze C : z= ^(w + w' 1 ) with w = pe w , 9 £ [0,2vr]|, (2.2) 

where C is the set of all complex numbers, and i = y—l is the complex unit. It has the foci 
at ±1, and the major and minor semi-axes are 

a=\(p + p- 1 ), b= 1 -{p-p' 1 ) 1 (2.3) 

respectively, so the sum of two semi-axes is p. The perimeter of £ p has the bound 

L{£ p ) < iry/p? + p-*, (2.4) 

which overestimates the perimeter by less than 12 percent (cf. |30j ). The distance from £ p 
to the interval [—1,1] is 

d p = \{p + p- 1 )-l. (2.5) 

We see that d p increases with respect to p, and d p —> + as p —} 1 + (so the ellipse reduces 
to the interval [—1,1]). Thus, by the theory of analytic continuation, we have that for any 
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analytic function u on [—1,1], there always exists a Bernstein ellipse £ p with p > 1 such 
that the continuation of u is analytic on and within £ p . Hereafter, we denote by 

A p := {u : u is analytic on and within £ p }, 1 < p < p max , (2-6) 

where £ Pmax labels the largest ellipse within which u is analytic. In particular, if p max = oo, 
u is an entire function. 

Throughout this paper, the Jacobi polynomials, denoted by J%'"(x) (with a, (3 > —1 and 
x £ I := (—1, 1)), are normalized as in Szcgo [33], i.e., 

1 J^{x)j:/{x)^{x) dx = 7 ^<W, (2-7) 



'-1 

where oo a ^{x) = (1 — x) a (l + x)@ , <J m , n is the Kronecker delta, and 

fi = 2"+Wr(n + a + l)r(n + /3 + l) 
7 " (2n + a + J 3 + l)n!r(n + a + /3 + l)' l ' J 
In Appendix [XJ we collect the relevant properties of Jacobi polynomials. 

In the analysis, we also use the following property of the Gamma function, derived from 
HJ Eq. (6.1.38)]: 

T(x + 1) = V2^x x+1/2 exp x + — J, Vx > 0, < 6 < 1. (2.9) 
Lemma 2.1. For any constants a, b 1 we have that for n>\,n + a> \ and n + b > 1, 

iXn+ry n ' (2 - 10) 

where 

T » ^ CXP l 2 (n + b-l) + 12(n + a-l) + ^^J- (2 ' U) 
Proof. Let 6*1,6*2 be two constants in (0, 1). We find from (|2.9p that 
Y{n + a) _{n + a-l) n+a - 1 ' 2 / 6>i 6* 2 



— — . exp - a + 6 + 



T{n + b) (n + b- \)n+b-i/2 *■ V 12(n + a-l) 12(n + 6-l) 

< (n + o-l)°- 6 (l+ 1 — 7 exp l-a+b+ 

\ n + b — 1 / V 



n + b-lJ V 12(n + a-l) 

^ a-bfi a-b\ a - b ( (a -b)(n + b- 1/2) 1 

<n a b H exp - a + & + v 



a - b , , ( a ~ 5 ) 2 A . = T o,6„, 



n ) \ n + b -I 12(n + a-l) 

< ?i a ~ b exp ( - 

F V2(n + 6-l) 12(n + o-l) 

where we used the fact that 1 + x < e x , for real x. □ 

Remark 2.1. Applying (|2.11j) to 7"^ leads to that for a, (3 > -1, n > 1 and n + a + [3 > 0, 

90+^+1 

2?1 + a + p + 1 

Note that for fixed a and b, 

T»' 6 = 1 + 0(0, (2-13) 
as it behaves like e 1 /™. □ 
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2.2. Main tools. 

Our starting point is the following important representation. 

Lemma 2.2. Let } be the Jacobi polynomial expansion coefficients given by 
1 rl 

In J -I 

If u £ A p with p > 1, we have the representation 



<' = —l u(x)J^(x)u a ^(x)dx, a,[3>-l, n>0. (2.14) 



-f I -^Ardz, n>0, (2.15) 
m J^> Jen wn+3+ 



where z = (w + w 1 )/2 with w = pe 10 , 9 £ [0, 2ir], and 

Ci = 4l t U n+j (x)J^(x)u;^(x)dx, n,j > 0. (2.16) 



a n,j a 8 

7n J -I 

Here, U n +j(x) is the Chebyshev polynomial of the second kind of degree n + j (cf. (|A.5[0 . 

Actually, the formula (|2.15[) - ([2.16[) can be obtained by assembling several formulas in 
Szego [33j, and then using the generating function of Uk(x) (cf. [1]). For the readers' 
reference, we sketch its derivation in Appendix [Bj 

The establishment of sharp bounds heavily relies on estimating . The following ex- 
plicit formulas follow from (|2.16|) and some properties of Jacobi polynomials listed in Ap- 
pendix [X] We remark that the formula (|2.19[) can be found in various books e.g., [T2l I29| ■ 
while the formula (|2.20p is due to Heine (see [H]). We also highlight that the formula (|2.2ip 
for the general Jacobi case seems new. 

Corollary 2.1. Letn>0. 

(i) For a = [3 > — 1 {ultraspherical/Gegenbauer polynomial^, 

<j = 0, for odd j. (2.17) 

(ii) For a = f3 = 1/2 (Chebyshev polynomial of the second kind), 

1/2,1/2 0r (n+l)\ 1/2,1/2 n f „ • ^ i fota\ 

" i> + 3/2) ; < = f° r ^ 1 - ( 2 - 18 ) 

(hi) For a = f3 = — 1/2 (Chebyshev polynomial), 



a 



(2^T(n+l) 
-1/2,-1/2 _ J ruTZTTo^ ' f° r even J 



= t r(n+l/2) ' J J ' (2.19) 

[0, for odd j. 

(iv) For a = f3 = (Legendre polynomial), 

, 2n+ir(l + l/2) T(n + l + l) 
°'° J 2 rTT+T) r(n + Z + 3/2)' for even j-M, 



a n,J 



0, /or odd j. 



In this paper, we do not distinguish between ultraspherical and Gegenbauer polynomials. 
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(v) For general a, ft > — 1 (Jacobi polynomial), 

at p _ 0F(2n + a + ft + l)r(n + a + ft + 1) 
~ 2r(n + a + l) 

7 (2-21) 
(-l) m r(2n + j + m + 2)T(n + m + a+ 1) 

•^-^ m!(j — m)!r(n + m + 3/2)r(2n + m + a + ft + 2) ' 

Proof, (i) . The property (|2.17p is a direct consequence of the parity of ultraspherical poly- 
nomials. 

(ii) . For a = ft = 1/2, we find from (|A.5|) and the orthogonality ([2TT ]> -(|2~g ]) that 

2~r^T7^' 5j ' ' 

where i^o is the Kronccker delta. Working out the constant leads to (|2.18[) . 

(iii) For a = ft = —1/2, if j = 21, we have 

a n,2l 1 = 1 -1/2,-1/2 n+ 2/ + 1 j 1 T ^+ 2l + 1 ^ J ' nl/2 ' 1/2 ( x )( 1 - 1/2 dx 

™ 2 ^ T „ (x) Jr7 i/2,-i/2 (x)(1 _ , 2) -i/ 2dx m 2ViT(n + 1) 



7n 1/2 '- 1/2 J-i r(n+l/2) 

which, together with (|2~T7)) . implies ([2~T9|l . 

(iv) For a = ft = 0, we derive from [TTJ Eq. (14)] that 

<ii = H>o J%°(x)U n+2l (x)dx=^- J%°(cos6)sm((n + 2l + l)6)d6 



In 

2n+ir(/ + l/2) r(n + ; + i) 

~~2 r(z + i) r(n + ; + 3/2) : 



Z > 0. 



This yields ((2~2Tj|) . 

(v) The formula ((2~2Tj) follows from a combination of (j2~g|) . (|A.4[) and (|A.5[) . □ 



With the aid of Lemma |2.2[ we can derive the following estimate, from which our sharp 
bounds are stemmed. 

Lemma 2.3. For any u € A p with p > 1, we have that for a, ft > —1 and n > 0, 

I -a f)| ^ M f I afli 1 I a fli 1 \^ i a R a,B\ 1 \ ,„ „„. 

I u « I < 3r(k,o I + -k,i I + 3LW+a ~ bJ' ( 2 - 22 ) 

1 V r j=o 1 

where M = max 26 £ p \u(z)\ and {cr"'j} are given by (|2.16p . 
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Proof. Since z = (w + w l )/2 e £ p with w £ C p (cf. (|2.i p -(|2.2 p ). we can rewrite 
(j2~T5)) as 



rri ^ " J 7 C u;"+jH 
i=o 

_^ oo 



dw 



2m ^ Jn 'i 

3=0 



2m 

+ 



T n,0 



C, 



w n+j+i 2m 



3=0 JC p 



u(z) 



dw 



dw 



1 

2rri' 



a,0 



1 oo 



a,/3 

27Ti ^ V " n ' J '+ 2 
J=0 



,,n+j+3 



(2.23) 



dw 



Hence, we arrive at 



C0\ < M 2n P I o,/3| M 2?r P I a,P I M 27Tp y-v I a,/ 
™ ' — 2tt / 9™+ 1 ' a,h ° ' ^tt "™+ 2 I "™' 1 ' 9-tt nn+3 2-^i r«J 



M i 



<*,P I 
7 n,0 



27T p 

M I a,/3| 
n+ l 1^,1 | 



2tt p 11 



ill 

Q n+2 



a,/9 





J+2 
1 



a.fl I 



This ends the proof. 



(2.24) 



□ 



Observe from the proof that we split the contour integral on £ p into two parts on C p , 
which actually allows us to take the advantage of cancelation of o~n'j+2 ~ a nj ■ Indeed, the 
bound (|2.22p is tight, as we will sec shortly that this argument can recover the best estimate 
for the Chebyshev case (see [3lJ Theorem 3.8] and (jl.lj) ). and improve the bounds in [38] 
(see {L2J)). 



2.3. Main results. For clarity of exposition, we first present the result on the general Jacobi 
polynomial expansions, followed by the refined results on the Chebyshev-type expansions 
(a = k — 1/2,(3 = I — 1/2 with k, I G N := {0, 1,2," •}), and Legcndre-type expansions 
(a = k,P = l with k,l E N). 

2.3.1. General Jacobi expansions (a, f3 > — 1). 

Theorem 2.1. For any u € A p (with p > l),a, /3 > —1 and n > 0, we have 

M 



\u a J\ 



< 



<*,0\ 
T n.O 



>n,l I 



'7o 



P P(P - 1) V 



where 



and 7™'P is defined in (J2T8J) . 

In particular, if a = ft, we have 



a ,p _ (2n + l)!r(7i + a + ;9 + l) 

r ™'° ~ 2 r(n + 3/2)r(2n + a + ^ + l) : 



a .p (/3-a)(2n + 2) a 
o~.„ i = — — —a, 



< 



M 



2n + a + ft + 2 



'To' 



71,0 ' 



(2.25) 



(2.26) 



1 V 7t7 



(2.27) 



8 X. ZHAO, L. WANG & Z. XIE 

Proof. By (12T221) . 



n I - - n \ a n,0 | + „+i 1^,1 | + „ +2 r»J+2 °n,j | 3 ■ (2.28) 



The factors <r"'g and <r"'f in (|2.26p are computed from (|2.2ip directly, so it suffices to 
estimate the infinite sum in (|2.28j) . Recall the identity (cf. [29]): 

U k (x)-U k _ 2 (x) =2T k (x), k>2. (2.29) 

Then we infer from (|2.f 6p that 

<f+2 " <f = f (Un +J+ 2(X) U n+j (x))j^(x)^(x) dx 

ln I (2.30) 
= 4fl/ T n+J+2 (x)J^(x)u J ^(x)dx, n,j>0. 

In J-l 

Thus, using the Cauchy-Schwartz inequality, the orthogonality (|2.7|) . and the fact |Tfc(x)| < 
1, leads to 



Kf +2 -<f|<-^=( / I^ i+a (*K''(*)<fa) ' <2W^. (2.31) 



Therefore, the bound (|2.25p follows from Y^Lo P J = 1/(1 — p I ),as/?>1. 
For a = /3, since |cr"'2/+il = 0? f° r all / > (cf. Corollary 12.11 (i)), we have 



El a . a. a , a I - \ ^ I a , a a , a I <r" O / ^0 

J=0 P i=0 P V 7 ™ 1 P 

This yields the refined bound in (j2~2T)) . □ 

Remark 2.2. Using Lemma 12. 1[ we can characterize the explicit dependence of the upper 
bounds in (|2.25[) and (|2.27j) on n, a, (3. Indeed, for a,(3 > — 1, n > 1 and n + a + (3 > 0, 

a g ¥I?m ^T(n + a + + l) (2n+l)! 



™' u 2 r(?i. + 3/2) T(2n + a + /3 + l) 

V^E a+/3+1 ' 3 / 2 n a+|8+1 ~ 3/ ' 2 ) (T2' a+,8+1 (2n) 2_ ( a!+/3+1 )) (2.32) 

= ^^W^ 1 ™ ^"(i + ofO), 

which implies 

„fs \2.2§\ \a~ /3\(2n + 2) a>j} \a - P\{2n + 2) WL3/2 2;tt+g+1 

1 "'^ 2n + a + /3 + 2 "'° " 2n + a + p + 2 2 a +? n 2n ^ ^ 



("-^^(l + OCn- 1 )). 
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Similarly, one verifies 

7^Jph 9 , , g , 7o^ n!r(rc + a + /?+!) 
7 «.0 ~ 1 +a + p+ j 2«+/ 3 +ir(n + a + l)r(n + /? + l) 

< (2n + a + /? + i) Mlff T WyW (2-34) 
1 (a + p + 2 J 

Consequently, we infer from the estimate (|2.25|) that for fixed a, ft > — 1 and n 1, 



ln " V 2Q+/3 ^ P ' V r (<* + + 2 ) P(P-1)J P n 

and likewise, we find from (|2.27[) that 

\uT\ < C nM (4 + 7 r(Q + 1) 2 f) (2.36) 
ln \2 2q v /r(2a + 2)p 2 -V P" V 7 

where C*„ = l + O^- 1 ). □ 

Remark 2.3. It is worthwhile to show that the bound obtained in this way is tighter than 
(fj~2j) obtained in [38]. Indeed, it follows from (|A.5|> . (|A.6bj) and (|2~7| that for n > 1 and 
i = 0, 1, 



In 

n+j 

T k (x)J^(x)u a ' (x)dx 



n+j -i r i 



2 

A:+?i+j + l odd 

= 4p f T n+] (x)J^(x)^(x)dx, 

In J -I 

where cq = 2 and Cfc = 1 for fc > 1. Following (|2.30l) - (|2.31l) . we have 



|<f|<2y / ^, n>l, j = 0,1. 



Finally, a straightforward calculation leads to 



M 



^| x K'il 2 7o ^ 



< 2M f 2 37^ 

p 1 p(p-i)V^ J -p"- 1 (p-i)V^" ' 

Moreover, we claim from (|2.27[) that the strict inequality holds, when a = > — 1. One 
may refer to Section U for numerical evidences. □ 

2.3.2. Chebyshev-type expansions (a = k - 1/2,0 = I — 1/2 with fc, Z € N). 

In view of (|2.19[) , it follows from (|2.23j) that the Chebyshev coefficient takes the simplest 
form: 



-1/2.-1/2 

f 1/2,-1/2 = I 3?L dw , (2.38) 

27ri J c w n+1 
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Thus, using (|2TT9"|) and leads to 



.,.,-1/2,-1/21 < V^T(n+l) M 
1 " 1 " r(n + l/2) 



(2.39) 



This leads to the estimate for the expansion coefficients, denoted by {u^} as before, in terms 
of {T n {x)} : 

2M 

K\<— , n>0, (2.40) 

as documented in e.g., |31j . 

For the second-kind Chebyshev case, we find from (|2.15|) the closed-form formula like 



1/2,1/2 



;l/2,l/2 _ 



u(z) 
„,n+l 



dz. 



but the contour integration is on £ p . It follows from (|2.18p and (|2.23|) that 



,1/2,1/21 < 



1 (n+l)l 



20Fl> + 3/2) 



,„n+l 



< 



(n + l)l A// 1 , 1 
2 r(n + 3/2)p" 



1 + — n 



(2.41) 



(2.42) 



Like (|2.40p . if we re-scale the expansion in terms of {U n }, i.e. 



u(x)U n (x)y/ 1 — x 2 dx, 



then we find from (|A.5|I and (f2~42]l that 



|fi?l = 



2 r(n + 3/2) 1/2 , 1/2 M/ 1 
^ r(n + 2) 1 " 1 S -» 1 + - 2 



P" 



(2.43) 



Remark 2.4. It is seen from (|2.41[) that the second-kind Chebyshev coefficient takes the 
simplest form on the contour £ p . This motivates us to estimate the contour integral directly 

by 



dz 



< 



P 



M f , , , M 



which implies 



By (ED), 



M L(£ p ) 

rM + l — 



(2.44) 



np 



< \ I 



1 1 

— < 1 + — ■ 

P 4 



Therefore, the estimate (|2.44[) is slightly sharper than (|2.43p . 



□ 



Some refined results can also be derived for a = fc+1/2,/3 = 1 + 1/2 with k, I £ N. Indeed , 
we find that {af;+ 1/2 -' +1/2 } can be computed explicitly by the following formula. 

Proposition 2.1. For any k,l,n,j G N, 

n+k+l 



fc+l/2,;+l/2 
7 ■ 

n,3 



1 \^ W fe + 1 /2,i + l/2 A /^l/2.1/2 A / 24 rx 

2 fc+l/2,! + l/2 d ™ VT™ <W+J> ^ Ab ) 

^ 7n m=ra 



where {d 



fc+l/2,J+l/2i n+fe+Z 



} = are ff™ 671 * n (|A.3[) . and 5 m) „_|_,- is i/ie Kronecker delta. 
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Proof. Using (|A.3|) (with a = (3 = 1/2), (|2.16[) and the properties of Jacobi polynomials (cf. 
([277)1 and (fOt ). leads to 



fc+l/2,i + l/2 



fc+1/2 , l+1/2 E C 1/2 ' z+1/2 / w*)# a,1/ ^)(i-^) 1/9 «k 



m— n 

-ri+fc+2 



fe+ 
7" 



This completes the proof. 



i — n-t-re+i , 

_t rL sr jhi/2,(+i/2,/ 1/2,1/2, 

1/2J+1/2U2 2-^ m V ' m °m,n+j, 



(2.46) 



□ 



Equipped with (|2.45p , we can obtain the bound for Chebyshev-type expansion coefficients 
by computing {d 1 ^ 1 ^ 2 ' 1 ^ 1 ^ 2 } explicitly. To fix the idea, we just consider the case: k = 1 
and I = 0. One finds 



and 



3/2,1/2 _ V7T n\(n + 2) 3/2.1/2 



.3/2,1/2 _ 1 ,3/2,1/2 _ 2n + 2 
' "+ 1 ~ 2n + 3' 

V^F (n + 1)! 



n -° 4 r(n + 3/2)' 

The estimate (|2.14p reduces to 

f 3/2,1/2 < ^£ 

n — n 

_ M 
- 7^7 

Thus, we have 



4 r(n + 3/2)' 



3/2,1/2 3/2,1/2 

3/2,1/2 , "Vi.l , a n,0 

V" + " 



^, 2 ' 1/2 = 0, 3 > 2. 



3/2,1/2 



1 



3/2,1/2 °n,l 



3/2,1/2-, 



y n,0 



|^/V/2| < 



/tF (n + 1)! M 



4 T(n + 3/2) 
and by (|2.10|) . we have for n > 0, 



lwn + 2 1 

! + ^r 7 + - 

/ \ n + 1 p 



(n+1)! 



< \fn exp 



8?i + 7 



1 



(2.47) 



(2.48) 



r(n + 3/2) - v " Vl2(2n+l)(n+l) 4n> 

Actually, the infinite sum in (|2.22l) does not appear for the Chebyshev-type expansions, 
which allows us to derive very tight bounds. However, for the Legendre-type expansions, 
some care has to be taken to handle this sum. 

2.3.3. Legendre-type expansions (a = k,[3 = 1 with k,l € N). 

We first consider the Legendre case. By (|2.8j) and (I2.2G|) . 



^0,0 

In 



2n + 1 ' j 
so the estimate (|2.27[) reduces to 



0.0 

7o 
0,0 



2n+l, a 



,o _ y / 7rr(n + 1) 
" r(n + l/2) ' 



-,0,0 1 



< 



M 
l>" 



V^r(n + 1) 2V2n+ 1 



r(n+l/2) 



1 



(2.49) 



In fact, we can improve this estimate, as highlighted in the following theorem, by using 



the explicit information of a. 



n,2V 



12 X. ZHAO, L. WANG & Z. XIE 

Theorem 2.2. Let be the Legendre expansion coefficients of any u <G A p with p > 1. 

Then for any n > 1, 



K°l^(^^Ml^> (-0) 

Proof. A straightforward calculation from (|2.20p yields 

o,o _ o,o n + 2? + 2 o,o , ^ n fo 

<r n ,2l+2 °n,2l- 2 (Z + l)(n + / + 3/2) f7 ™< 2 " " ' ^ 5iJ 

which implies {c^;} is strictly descending with respect to I. Hence, we have 

o,o ^ n + 2 o o 

^2/ ^ 7T— T CT «.0> ( 2 - 52 ) 



0,0 _ 0,0 I ?z + 2Z + 2 o,o < n + 2 ,o 

J n,2l+2 °n,2l\- 2 (l + 1) („ + / + 3/2) ^ " 2/1 + 3^°' 



where we used the fact that n + 21 + 2/((l + l)(n + i + 3/2)) is strictly descending with 
respect to I. Then, we obtain the improved bound from (|2.22l) : 

I -ooi M 0,0/ n + 2 i \ 0rM r(n+l) / n + 2 1 



<-<:o i 



p" ™< U V 2n + 3 p 2 <+ 2 / p" r(n+l/2)V 2n + 3 p 2 - 1 



and by (|2T0|) . 

T(n+1) _ / 8n-l \ .„ 

i^+T72T^ cx KiM2^)' n ~ ^ 

This completes the proof. □ 

Remark 2.5. We compare the bound in (|2.50p with the existing ones. Davis [12] Page 313] 
stated the bound 

i-ooi < 2n+l ML(£ P ) IO 2n + 1 ii\J p 2 + p~ 2 M 

where clearly the algebraic order of n in the numerator is not optimal. The following 
asymptotic bound can be obtained from [27j Eq. (32) and Eq. (38)] and QU Eq. (12.4.25)]: 



.0 01 / My^m vV + 1 

K'° < „ 2 1 ' n>1 ' 

p™ p — 1 



while the asymptotic estimate derived from (|2.50l) is 

,o,oi ^ M^ip 2 - 1/ 2 
p" p 2 - 1 

which is sharper. Another bound for comparison is obtained in the recent paper 



o. 0| VnM/ _L_ 
™ 1 S p™ V i+ p 2 -iy 



n>l, (2.56) 



which is also inferior to our estimate (|2.50[) . Some comparisons in numerical perspective are 
given in Section [U □ 

Like the Chebysheve case, we can derive similar refined estimates for Legendre-type 
expansions with a = k, (3 = I and k, I € N. The counterpart of Proposition 12.11 is stated as 
follows, which can be obtained by using (|A.3|) (with a = f3 = 0), (|2.16|) and the properties 
of Jacobi polynomials (e.g., (|2.7[) ) as before. 
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Proposition 2.2. For any k,l,n,j £ N, 



n+k+l 

k,l 1 ,k,l 0,0 0,0 c 7 \ 

T n m—n 



where {^m'}^-n~ are the same as in ()A.3[) . and {ff^ U n+:) _ m } are computed by (|2.20[) . 

Once again, to fix the idea, we just consider the case: k = 1 and I = 0. One finds 
d^ a = l,d„' +1 = — 1, and 

1,0 / o.o o.o _ 0,0 o.o \ _ n + 1 ,o _ n + 1 ,o 

°n,j — 1,0 l'» a n,j 7n+l a n+lj-l) ~ 2n _|_ 1 °™ J 2 n + 3°' n + 1 ^- 1 " 

By Q, 

1,0 _ n + 1 o.o 1,0 _ n + 1 o,o / > n 
"• 2i ~ 2n + 1 CT "< 2/ ' CT "' 2i+1 ~ 2n + 3 CT ' l+1 - 2 " ~ 

Therefore, with (f2~5TT) and (|2"32|) . the estimate (|2T2"2]) reduces to 

I -1,01 ^ ^ \ 1 ,0 1 |^ | 1,0 I , \ ~* I 1.0 _ 1.0 1 1 

\ U n I - „ \ a n,0\ + n+1 \ a n,l \ + „ +2 l CT n,J+2 I 

F H H j=0 H 

_ M n + 1 / o.o 1 | o.o _ o,o I J_\ 

" p« 2n + 1 V ™'° + p2 2_, 1^,2^2 CT n ,2/| p 2i J 

^ n + 1 / o,o J_V^| o.o o.o ij_ 

2n + 3 v p 2 ^ l a "+ 1 > 2 i+ 2 (T «+i,2il p2Z 

oo n+1 M / 71 + 2 1 \ 00 n + 1 M / n + 3 1 
< cr ' 



1 M / 

'^n+l^l 1 ^ 2?i + 3p 2 -iy ^"" +1 >°2n + 3p»+ 1 ^ 

.0 

+ 

^' I < 



Working out the expressions of cr^'o and o~ n ', i by (|2.26|) . we have 



M ^T(n + 2) fl n + 2 1 1 n+1 / n + 3 1 \\ 

p n T(n + 3/2) I 2 2(2n + 3) p 2 - 1 p 2n + 3 V 2n + 5p 2 -l/i' [ ' ' 

Note that the ratio of the Gamma functions can be bounded as in (|2.48l) . 

The same process applies to other k, I £ N, but the derivation seems tedious. 

2.4. Estimates for truncated Jacobi expansions. Given a cut-off number jV > 1 and 
N £ N, we define the partial sum 

JV-l 

{-K a /u){x)=Y J K' J^^), (2-59) 

n=0 

where {u"'^} are the Jacobi expansion coefficients defined in (|2.14|) . To this end, let L 2 aiP (I) 
be the weighted L 2 -space on I = (—1,1), and its norm is denoted by || • H^q./s, where we 
drop the weight function, if a = (3 = 0. 

Notice that -n^u is the L 2 Q j)9 -projection of u upon P/v-i (denoting the set of all algebraic 
polynomials of degree at most N — 1), that is, tt^ u is the best approximation to u in the 
norm || • || w a,/3. With the previous bounds for the expansion coefficients, we can estimate the 
truncation error straightforwardly. 
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Theorem 2.3. For any u € A p with p > 1, and a, (3 > — 1, we have 



a,0 I 
TT N U — U\ 



< 



\a-0\ 



2V^ 



C N M 



,iV-l 



2 a +? V ' p ) ' p(p- 1) 
where 7q is given in (|2.8p and Cjv « 1. 

Proof. By the orthogonality (cf (|2.7[) - (|2.8[) ) of Jacobi polynomials, we have 



(2.60) 



7Tjy U — M 



It follows from the estimate of \u%'P\ in Theorem 12. 1[ and a combination of (|2.12[) - (|2.13[) 
and ([2~52>([2331 that for n > AT > 1, 

2 



,<*,/? < 



< 



M 


1 




r 


p n 




C n M 


p n 



a, 8 



a, 8 



1| a,B\ I a, 8 



where C„ = 1 + 0(n L ). Therefore, we have 



ct,/3 



7Tjy It — It 



'n,l 



P(P-1) 



7o 



2^ 



7o 



a,/3 



< 



2"+^ v p 



p(p-l) 



n=N H 



-JV-l 



This ends the proof. 



□ 



Remark 2.6. Note that {4-[ J" ,/3 }n>; are mutually orthogonal with respect to u a + L P+ l ; so 



we can estimate (71^ u — u) 



in a similar fashion. 



□ 



+1,0+1 

Remark 2.7. Some refined estimates can be obtained from the refined bounds for special 
cases, e.g., oj = /3ora = /3 = 0, —1/2. Here, we just state the result for the Legendre case: 

1 \ CVv^FM 



I o.o 1 1 r 

1 TT^ U - U|| < I 1 



2(p 2 -l)V pN-i^TZTi 



(2.61) 



where Cn ~ 1 as before. It follows from Theorem 12.21 and the above process. Note that 
Xiang [35] derived the following estimate for the Legendre expansion: 

2y/2M 



0,0 II ^ 

ir N u-u\\ < ^ 



p"- 2 (f-l 

The estimate (|2.6ip seems tighter than this one. 



(2.62) 

□ 



3. Error estimates for Gegenbauer-Gauss quadrature 

3.1. Preliminaries. The Gegenbauer-Gauss quadrature remainder (|1.3[) - (|1.4[) with the nodes 
being zeros of the Gegenbauer polynomial J" ,a (x), takes the form 



E G n G [u] 



hi 



£ D Jn (Z) 



u(z) dz, Vm G A p , 



(3.1) 
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where Q"' a (z) is defined as in (|B.2[) . namely, 

w)- J ^fs3&p**'w±^<>v±^. ( 3. 2 , 

J ~ 1 ]=Q 1=0 

As already mentioned, the analysis of quadrature errors (even for the Chebyshev case) 
has attracted much attention (see e.g., [UJ [TO] EJ [TT1 US] [16j 123 [26] ) . Just to mention that 
Chawla and Jain [TTJ Theorem 5] obtained the estimate: 

l B » G M|^^rrp v,ei„ vn>i, (3.3) 

Hunter [35] derived the general bound 

K G M1< ^n-2 (p 2_ 1) . n>l, (3.4) 

and some refined results for a = ±1/2 and /? = ±1/2 by expanding <5"' a /J"' Q into the 
Laurent series of w in the disk enclosed by C p , and manipulating the series. It is worthwhile 
to note that Gautschi and Varga [17] estimated the Jacobi-Gauss quadrature (with J™'" 
and Q"'* 3 in place of J"' a and Q"'" in (|3.1[) . respectively) by 

|2tf G M| < Tr-S^ML^maxlQ^^/J^^)!, (3.5) 

and attempted to find the exact maximum value on the Bernstein ellipse, which was feasible 
for a = ±1/2 and = ±1/2 again. Some conjectures and empirical results were explored 
in [17] for the general Jacobi case. 

Using the explicit expression of Legendrc polynomials on the Bernstein ellipse (see e.g., 
[T2l Lemma 12.4.1]), Kambo [27] obtained the bound for the Legendre- Gauss quadrature: 

tpLG r 1 1 -1 0.0,,r/c , ™*ze£ P \Qn°( z )\ ^ d n M p 2 + 1 /- 

mm ze£p | J„ (z)| P P - * 

where < d„ < n. While this bound is only valid for p > it holds for all n, when 
compared with the asymptotic estimate (with u» 1) for the Legendre-Gauss quadrature 
in [9]. 

In what follows, we aim to extend our analysis to estimate E GG [u] in (|3.1I) . The essential 
tools include the explicit formula for the Gegenbauer polynomial J"' a (z) on £ p derived in 
our recent paper [39], and the previous argument for estimating Q^ a (z). Let us recall the 
important formula stated in (391 Lemma 3.1]. 

Lemma 3.1. Let z = ^(w + w^ 1 ). Then we have 

n 

JZ' a {z)=A a n Y,9t9n-k^ n - 2 \ n>0, o>-l, a? -1/2, (3.7) 

k=0 

where 

y ° JL fc!F(a + l/2) ' - - ' n r(a + l)r(n + 2a+l) V ; 

Remark 3.1. This formula excludes the Chebyshev case. For a — —1/2, we define 

So" V2 = -9„- 1/2 = 1. V2 = 0. 1 < * < n ~ 1, and = T ^t}{ 2) , (3.9) 
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since (see e.g., [T2]) 



T n (z) = -(w n + w- n ) 



2A, 



: r-i/2,-i/2/„i 

-1/2 ™ 



(3.10) 



Hence, we understand that (|3.7[) holds for a = — 1/2 with the constants given by (|3.9p . □ 

3.2. Main results. We adopt two approaches to estimate the quadrature remainder. The 
first one is to expand <3"' a /</"'" m Laurent series of w € C p , and then we use an argument 
as for Theorem 12.31 to obtain the tight error bound. However, this situation is reminiscent 
to that in Gautschi and Varga |17j . that is, computable bounds can be derived for general 
a. We highlight that the computational part (see p.lip ) is independent of p and u. 

The second approach is based on an important relation between the quadrature remainder 
and Gegenbauer expansion coefficient (see p.22p ). 

The main estimate resulted from the first approach is stated as follows. 



Theorem 3.1. For any u € A p with p > 1, we have that for a > — 1 and n > 1, 



E 



GG 



«] < t£ 



1 



I A*n,0 | + | M„,2i+2 M n ,2J | p2 _l 

where {/U„'2;}i>o ar c computed by the recursive formula: 



M 

r,2n ' 



(3.11) 



a, a min{n,/} 
a,a _ _J_ f a n,2l 



4 a 



E a ot a. a \ ; \ i a. a 



fe=i 



"n,0 
A a n a 



(3.12) 



Proof. A straightforward calculation from (|3.2|) (note: cr"'^ , j = for all I > 0) and p.7p 
leads to 



EOT,, 2/ a, a a a \ A a a a, a 

,.,2n+2i+l Wlth a n,2l = A n 2^ 9k9n-kVn,21 



min{n,Z} 



r 

■J a 



2l~2k ' 



(3.13) 



A;=0 



so solving out /x™'^ yields p.!2p . 

Next, following the same lines as the derivation of (|2.23[) . we infer from (|3.ip and p,13p 
that 



2?rp 



1=0 

a, a a, a I 

" 2tt L 2n+l l^n.O I + 2n+3 l^n,2i+2 ~~ ^n,2l 1 

H H 1=0 



1 1 

n2l 



'« p2n 



Mn.O 



1 X I a, a 

3 l^n,2i+2 
H 1=0 



a, a 
^n,2l I ^27 



M 



n p 2n 
p2n 



a, a 


+ max 

i>0 




a 

21+2 


a.Q: 
- M„,2/ 


a.a 
Mn,0 


+ max 

i>0 


a 


a 

2Z+2 


a.Q 
" Ai„,2/ 



„2;+2 



1=0 



(3.14) 



This completes the proof. 



□ 



JACOBI EXPANSIONS AND GEGENBAUER-GAUSS QUADRATURE OF ANALYTIC FUNCTIONS 17 

Remark 3.2. We find from ([3~T2"j) that for a = -1/2, 

-1/2,-1/2 _ 2tt I -1/2,-1/2 _ -1/2,-1/21 27t£ k ,q _ i/i.i \ 

Mn.O — -1/2-1/2' l^n,2/+2 I - -1/2-1/2' ^ ' HlOQ^t -f" 1,7&J, 

where <5 Ki q is the Kronecker delta. Hence, it follows from (|3.14p that 



2ttM 

E n G M\<^ + 7 E 7 U^j\=^- V ^>1, (3.15) 



1 1 



2wM 



which is the same as Q3.3P derived in [11] . □ 
Remark 3.3. We find from (f3TT2|) that for a = 1/2, 

1/2,1/2 J(- 1 ) K f 1/2,1/2 ' if «:=mod(l,n + l) = 0,l, 
<2i = S ^ 7n 7 (3-16) 

[0, otherwise, 

which implies 



oo oo 
V^l ,1/2,1/2 _ 1/2,1/2 1 1 I 1/2,1/2 1/2,1/2 | V^| 1/2,1/2 

/ y |M n ,2i+2 ^n,2; I 2i — l^",2 H'nfi \ "T I" 

(=0 ^ j=l 



2j"("+l) I p2j(ra+l)-2 



E/l 1/2,1/2 1/2,1/2 I 1 I 

Vl^n,2j(n+l)+2 _ M n ,2j(n+1) I 2j(n+l) + 1^ 



1/2,1/2 



)lp2j(n+l) ^ l^n,2j(n+l)+2 I /3 2j(«+l)+2 



„-l\2 



K 1 /o / -1\2V"^ 1 \ 7T 1 

9 1/2,1/2 l 2+ (P + P ) „2i(n+l) j = 9 1/2,1 



9 , (P±P_X 

2 1/2.1/2 V" ' VP ' ^ ; ^ fl2j(n+l) y 2 1/2.1/2 V 2 n+2 _ J 

In j=l r In ' 

Hence, it follows from (|3.14j) that for the Chebyshev- Gauss quadrature of the second kind, 



2p 2n y p 1 V pihH-2 _ 1 J J 2(p 2n + 2 - 1) 

Note that Hunter [551 (4-8)] obtained the following estimate by a delicate technique: 

\ E n M\ < 2(p 2 n +2 _ ^ • ( 3 ' 18 ) 

We see that (|3.17p is sharper. □ 

For general a > ~ 1, the derivation of an explicit bound for 

e^maxC,,, K,l ~ 7n' a K;2;+ 2 ~ infill n >^ ( 3 - 19 ) 

seems nontrivial. We have only empirical results bases on computation. Some indications 
are listed as follows. 

1 i 

(i) Observe from f|3 . 16[) that for fixed n, {9 ', };>o are (n + l)-periodic (see Figure [3~T1 

(a) ), and the maximum is attained at I — j(n + 1), j — 0, 1, • ■ • . We compute ample 
samples of n, I and a, and find very similar "periodic" behaviors (see Figure 13.11 

(b) -(c) for a = 0, 1). 

(ii) Another interesting empirical observation is that for fixed a, the maximum value 
9" converges to a constant value, and it decreases as a increases (see Figure 13.11 
(d)). Note that for the Legendre case, 0° w 4. 

Now, we turn to the second approach. The main result is summarized below. 
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n 



(d) 0" with various a 



FIGURE 3.1. (a)-(c): Profiles of 0™ , with n = 36, a = 1/2, 0, 1 and < / < 250. 
We mark by "□" the location of the maximum value O" is attached, (d): The 
maximum value O" with a = 0, 3/2, 5, 10 and 10 < n < 100, where we compute 
{0^} for I up to 1000. 



Theorem 3.2. For any u € A p with p > 1, and /or a > —1 and a ^ —1/2, rae have 

| F GGr 1 i < C n My/Z (yfF r(a + l) 2^2 \ / (1 + p' 2 )^ 1 / 2 , a > -1/2, 

1 " [ J| " P 2 " V2 2 « + v /r(2a + 2) P 2 ^ I/ \ (1 - p- 2 Y^/\ a<-l/2, ^ 

and in particular, for the Legendre case, 



Ei: G \u]\ < 



where the constant C n ~ 1- 
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Proof. We carry out the proof by using the important relation, due to (|3.1[) and (|B.1I) : 

(3.22) 



K G M| < - — 7 v. a , - 4 

1 ram ze s p \Jn [z)\ tt 



min zefp |J"' Q (z)| ' 

Since the numerator has been estimated in Theorem 12. II falso see (|2.36|) ). it suffices to deal 
with the denominator. 
By [391 (4.6)], we have 

\j^(z)\>c\A«\ n "" 1/yi fa+p- 2 )-"- 172 * if «>-V2, 

1 jl - n|A,l |r(a + 1/2)1 1(1- — 1/2, if a < -1/2, 

(3.23) 

2 2 >" [(l + p- 2 )-"- 1 / 2 , if a > -1/2, 
" " p- 2 )-*- 1 ' 2 , if a < -1/2, 

where C„ « 1. Note that in the last step, we dealt with |^4"| as 

IP |r(2a + l)| r(n + Q + l) = \T(a + 1/2)| T(n + a + l) > |r(a + l/2)| 
1 r(a+l) T(n + 2a + l) 2" 2q 0F r(n + 2a + l)~ " 2- 2a v ^ji a ' 

where we used Lemma 12. II and the property of Gamma function (see [1]): 

T(z)T(z + 1/2) = 2 1 - 2z V^r(2z). 

Therefore, a combination of (|2.12[) . (|2.36p and (|3.22[) - (|3.23[) leads to the desired result. 
Using the refined estimate (|2.50[) . yields (|3.21[) . □ 

4. Numerical results and comparisons 

In this section, we present various numerical results to show the tightness of the bounds 
derived in this paper, and to compare them with other existing ones mentioned in the 
previous part. 

In the first example, we purposely choose the Chebyshev and Lcgcndrc expansions with 
known expansion coefficients: 

n— 1 v n— 

which follow from generating functions of Chebyshev and Legendre polynomials (cf. [33] ). 

Note that the function Ui has a simple pole at z = 5/4, so the semi-major axis (cf. (|2.3[1 ) 
should satisfy 

1< a = (p + p- 1 )/2 < 5/4 =>■ l<p<2. 

One also verifies that 

M = T£ lui{z)l= (2 P -i)(2- P y 

Then the estimate (|2.40j) reduces to 

Un = 77— T < 777 77777 s T-= B n(p), 1 < P < 2 , Tl > 1. 

2 „-i - (2p-l)(2-p)p™- 1 ™ ^ r 
Similarly, for the Legendre expansion of U2, the result (|2.50|) becomes 



1 ^ y/nn ( ^ t n + 2 1 ^ ^ 8n — 1 ^ / 4,0 



V - F - 7^ V 1 + 2n + 3 p 2 - J ° XP i 12n(2n - 1) J V (2p - 1)(2 - p) '= B " {p) ' 
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for 1 < p < 2 and n > 1. 

We take p = 1.98, and plot the exact coefficients and and the bounds 2?^ and 
in Figure 14.11 (a) and (b) , respectively. Actually, the bound for the Chebyshev case (see 
(jl.ip ) can be considered as one benchmark for illustrating tightness of the upper bound. 
Indeed, the result for the Legendre case stated in Theorem 12.21 seems as sharp as that for 
the Chebyshev case. 




n n 

(a) Chebyshev case (b) Legendre case 



FIGURE 4.1. Expansion coefficients of wi,U2 in (|4.1|l against their error bounds. 
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(a) e n (p) (b) e n (p) with p near 1 

FIGURE 4.2. (a): Comparison of error bounds for Legendre expansions in (|2,50[) 
and (|2.56f) . (b): Samples of e n (p) for p close to 1. 

Next, we compare the bounds for the Legendre expansion coefficients in Theorem l2.2l and 
(|2.56p (obtained by [38]). For clarity, we drop the common part M^fnj p n , and denote the 
remaining factors in the upper bounds by 
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In Figure B31 fa), we plot the difference e n (p) := b(p) — b n {p) for various p and 1 < n < 80. 
We see that e n (p) > and the difference is of magnitude around 6, when p is close to 1. 
Moreover, for fixed p, the difference is roughly a constant for slightly large n. In Figure B~2l 
(b), we plot some sample e n (p) for p close to 1, we see that our bound is much sharper. 




(a) a = l,/8 = (b) a = = 2 



FIGURE 4.3. (a): Comparison of error bounds for Jacobi expansion with a — 
1,P = in (|1.2p and (|2.58|) . (b): Comparison of error bounds for Gegenbauer 
expansion with a = /3 = 2 in (|1.2|) and (|2.27|) . 

We next make a similar comparison of bounds for Jacobi and Gegenbauer expansions. 
For example, for a = 1 and /3 = 0, we extract the factors in (|1.2p and (|2 . 58[) by dropping 
My/n/p n . Wc plot in Figure l4~3l fa) the difference of two remaining parts (i.e., that of (|1.2j) 
subtracts that of (|2.58|l ). Once again, our bound is much tighter. Likewise, we depict in 
Figure l4~3l (b) the extracted bounds from (|1.2|) and (|2.27|) with a = /3 = 2. The situation is 
mimic to the Legendre case, where the bounds obtained in this paper are sharper. 

Finally, we turn to the comparison of error bounds for the Gegcnbauer-Gauss quadrature 
remainder. For a = 1/2, we extract the factors in (|3.17j) and (|3.18j) by dropping M/p 2n 
as before. We plot in Figure 14.41 (a) the difference of two remaining parts in (|3.18|) and in 
(|3.17[) ). Once again, our bound is much tighter. Likewise, we depict in Figure B~4l (b) the 
extracted bounds from (|3.4[) and p. lip with a = 2, and observe similar behaviors. 

Concluding remarks 

In this paper, we derived various new and sharp error bounds for Jacobi polynomial 
expansions and Gegenbauer-Gauss quadrature of analytic functions with analyticity charac- 
terized by the Bernstein ellipse. Wc adopted an argument that could recover the best known 
bounds, and attempted to make the dependence of the estimates on the parameters explic- 
itly. Both analytic estimates and numerical comparisons with available ones demonstrated 
the sharpness of the error bounds. 

Appendix A. Jacobi polynomials 

Wc collect some properties of Jacobi polynomials used in the paper. For a, ft > —1, the 
Jacobi polynomials (see e.g., [33]), denoted by J"'"(x),x € I := (—1, 1), arc defined by the 
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(a) a = 1/2 (b) a = 2 



FIGURE 4.4. (a): Comparison of error bounds for the Gegenbauer-Gauss quad- 
rature with q = 1/2 in (|3.17p and ()3.18|) . (b): Comparison of error bounds for 
the Gegenbauer-Gauss quadrature with a = 2 in (|3.4f) and (|3.1ip . 



Rodrigues' formula 



The Jacobi polynomials satisfy 



(l-x)JZ +1 >f>{x) = 



(-1)™ d n 



2 n n\ dx r 



(l-x) a+n (l+x) 



P+n 



n > 0. 



2n + a + (3 + 2 



(l + x)J^ +1 (x) = 



2n + a + f3 + 2 

As a direct consequence of (|A.2[) . we have that for any fe,(sN = {0, !,•••}, 



n+k+l 



(1 - x) k (l + x) l J% +k >P +l (x) = Y, d^ +w J^(x), 



(A.l) 



((n + a + l)J^{x) (n + 1) J^x)) , (A.2a) 
; ((n + p + l)J^(x) + (n + l)J a n Ux)). (A.2b) 



(A.3) 



where {'i" +fe '^ + '}™^ r f +i is a unique set of constants (with d"'^ = 1), computed from (|A.2[) 
recursively. Here, we sketch the proof of (|A.3[) . To this end, let {cj} be a set of generic 
constants. Using (|A.2a[) and (|A.2b|) repeatedly leads to 

{l-x) k (l + x) l J% +k 'P +l (x) 

= (i - x f-\i + x) \ Cl (x) + C2 j^- 1 ^)) 

= ■ • • = (1 + X) 1 Y CmJ^ +l (x) =■ = Y C >nJ^(x). 



This yields (jA~3)) . We point out that for a = /3 = 0, {(1 - x) k (l + x) l J^ 1 } (up to a certain 
constant factor) are defined as generalized Jacobi polynomials in [2"2"| . 
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The following formula, derived from ^ Lemma 7.1.1] (also see [3H Theorem 3.21]), was 
used for the derivation of (|2.21l) : 



®(a,l3,a,b) 



1 



1 



J a 4 j {x)J^{x)^{x)dx 



T(n+j + a + l) (2n + a + /3 + l)r(n + q + ^ + l) 
T(n + j+a + b+l) T(n + a + 1) 

^ (-l)™r(2n + j + m + a + b+ l)T(n + m + a + 1) 
^-^ m\(j — m)!r(n + m + a + l)T(2n + m + a + (3 + 2)' 



(A.4) 



for a,b,a,0 > —1 and n, j > 0. 

Let T n (x) = cos(n arccos(x)) be the Chebyshev polynomial of the first kind of degree n. 
Then the second-kine Chebyshev polynomial, denoted by U n (x), can be expressed by 



U n {x) = 



3in ((n + 1) arccos(x)) T' n+l (x) 



n + l 



2 / 1/2,1/2 
In 



The Chebyshev polynomials enjoy the following important properties: 

J-^-^(x) = j-^'\i)T n {x) = r( " + Y 2) r K (x), 



n-1 



1 

T' n {x)=2n J2 ~W 



k=0 
fc+n odd 



where cq — 2 and = 1 for k > 1. 

Appendix B. Proof of Lemma |2~21 



We first show that 



1 



where 



1 



and 7^ ,/3 is given by 



2-) 



a,/3 



1 ^^(iK'^l) 



z — a; 

Recall the Cauchy's integral formula: 

u(z) 



d n , . n! 
- — u(x) = — — 



dx r 



2ni 



(z - x) n+1 



dz. 



Using the Rodrigues' formula (|A. 1|) and integration by parts leads to 



7- 



ct,/3 



u(a;)J n a ' ? (s)w a ^(a;) dx 



-l 



(-1)" {\«+n, 0+ n {x) d^ u{x)dx 

J 



1 



a,/3 2 n n\ 



E3J i i 



dx r > 
u(z) 



7 «,/3 2 «n! 7_! V27ri J £p (z-x) 

(x) 



'')—dz\uj a+n - fs+n {x)dx 

\n+l ) \ ' 



1 1 



1 , ,a+n,/3+n 



W-i V *)"+* dx ) u ^ dz 



(A.5) 

(A.6a) 
(A.6b) 



(B.l) 
(B.2) 

(B.3) 



(B.4) 
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We find from integration by parts that 



1 a+n,j3+n( \ (—T\ n f 1 1 rl n 

Til ] . v — f d'r-n * i ^ i 



_i (z — a;) n+1 n\ J_ 1 z — x dx 

Inserting (|B.5j) into ()B.4j) . we derive from the Rodrigues' formula (|A.lj) that 



where <3" ,,3 (2:) is given in (|B.2j) . 

Since z = (u>+u; -1 )/2, we have from the generating function of the Chebyshev polynomial 
of the second- kind (cf. [T]) that 

12 1 2 ^ Uk ( x ) 



E^P- (B.6) 

z — ^ ill n- 



Z — X w w — 2xw 1 + 1 It) ^— ' ID 

fc=0 

Inserting it into (|B.2[) . we find from the orthogonality of the Jacobi polynomials (cf. (|2.7p ) 
that 



T™ fe=o 



oo „i 

" fe=o 

= ^E J+T f_Uk{x)J^{x)^{x)dx (b.7) 

7n „•_ n •'— 1 „•_ n 



where we defined 

<f = -is f U n+j (x)J^(x)^(x)dx, 
Substituting the last identity of (|B.7|) into (|B.1|) leads to the desired formula (|2.15[) . 
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